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® Platform Stability Analysis - aircraft attitude variations 


• Microburst Detection Without Conventional Ground Clutter Rejection 

-autoregressive modelling 
-microburst tracking 
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- adaptive noise cancelling 

- simulated microburst in real clutter data 
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Abstract 

Radar data collected during the 1991 NASA flight tests have been selectively 
analyzed to support research directed at developing both improved as well as new 
algorithms for detecting hazardous low-altitude windshear. Analysis of aircraft 
attitude data from several flights indicated that platform stability bandwidths 
were small compared to the data rate bandwidths which should support an 
assumption that radar returns can be treated as short time stationary. Various 
approaches at detection of weather returns in the presence of ground clutter 
are being investigated. Non-conventional clutter rejection through spectrum 
mode tracking and classification algorithms is a subject of continuing research. 
Based upon autoregressive modelling of the radar return time sequence this 
approach may offer an alternative to overcome errors in conventional pulse-pair 
estimates. Adaptive filtering is being evaluated as a means of rejecting clutter with 
emphasis on low signal-to-clutter ratio situations, particularly in the presence of 
discrete clutter interference. An analysis of out-of-range clutter returns is included 
to illustrate effects of ground clutter interference due to range aliasing for aircraft on 
final approach. Data are presented to indicate how aircraft groundspeed might be 
corrected from the radar data as well as point to an observed problem of 
groundspeed estimate bias variation with radar antenna scan angle. A description 
of how recorded clutter return data are mixed with simulated weather returns is 
included. This enables the researcher to run controlled experiments to test 
signal processing algorithms. In the summary research efforts involving 
improved modelling of radar ground clutter returns and a bayesian approach 
at hazard factor estimation are mentioned. 
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NOTES 


Roll angle variation with time during approach to 

L U i2 Way K 2 7 at PH ^' Data were record ed from DATAC 
wifteach frame of radar data at a frame rate of 29 25 
frames per second. 
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NOTES 


Frequency spectrum of roll angle time variation 
during approach to runway 27 at PHL. 

Data were recorded from DATAC with each 
frame of radar data at a frame rate of 29.25 
frames per second. 
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Crab Angle vs. Time PHL f3r3s2.m6 
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NOTES 


Crab angle variation with time during approach to 
runway 27 at PHL. Data were recorded from DATAC 
wifteach frame of radar data at a frame rate of 29 25 
frames per second. 
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NOTES 


Frequency spectrum of crab angle time variation 
during approach to runway 27 at PHL. 

Data were recorded from DATAC with each 
frame of radar data at a frame rate of 29.25 
frames per second. Crab angle mean was removed 
prior to spectral analysis. 
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Pitch Angle vs. Time PHL f3r3s2.m6 
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NOTES 


Pitch angle variation with time during approach to 
runway 27 at PHL. Data were recorded from DATAC 
with each frame of radar data at a frame rate of 29.25 
frames per second. 
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NOTES 


Frequency spectrum of pitch angle time variation 
during approach to runway 27 at PHL. 

Data were recorded from DATAC with each 
frame of radar data at a frame rate of 29.25 
frames per second. The mean value of pitch angle 
was removed before spectral analysis. 
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Orlando Flights Though Events 142,143,144 
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NOTES 

Ground tracks of three legs of the NASA flight through 
the microburst in Orlando on day 171. This event 
was numbered 141 on the first pass, 143 on the second, 
and 1 44 on the third. The aircraft was at about 1100 feet 
traveling in excess of 200 knots. The indicated position 
of the microburst is an estimate of the core position 
based upon radar data. The TDWR is at 0,0 on this plot. 
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NOTES 

During the flight through event 143 at Orlando a 
snapshot of the radar return looking into the 
microburst is analyzed using a Fourier transform 
of the I & Q sequence (96 samples) taken at an 
antenna azimuth of -0.25 degrees. The aircraft 
was on the track labeled 143 on the previous slide 
and located at about 8 km west and 22 km north of 
the TDWR. The range cell at the zero crossing of 
the "s-curve" characteristic is about range cell 35 
which is appoximately 5 km ahead of the aircraft. 

In the presentation the Doppler power spectrum 
has been thresholded and then a point density 
plot is used to indicate spectrum intensity. Doppler 
windspeed is on the abscissa. Range cells from 6 
to 96 are indicated on the ordinate. No clutter 
rejection filtering has been used. 


597 



Autoregressive (AR) Modelling of Radar Return 
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NOTES 

Linear modelling of the I & Q sequences as if they 
were the output of a linear all-pole model driven by 
white noise is being used to investigate the feasibility 
of detecting modes in the return Doppler spectrum 
without the use of clutter rejection filtering. In 
situations where the clutter is particularly strong 
or may tend to bias spectrum mean estimates even 
when attempts at clutter rejection are made, this method 
is viewed as a possible alternative. It also provides 
a method for estimating a spectrum as an alternative 
to the FFT. A second order AR model is comparable 
to the pulse pair algorithm in terms of processing load 
and can give a useful spectrum estimate for further 
processing. 








AR Based Velocity Estimates for Event 143 
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NOTES 


The snapshot of event 143 in Orlando at a time very 
near to that shown in the earlier scatter point plot is 
shown here after autoregressive modelling. Again no 
clutter rejection filtering has been used. The bubble 
center locations indicate the spectrum mode Doppler 
velocity estimates for each range cell. The size of the 
bubbles indicate relative mode strengftjjs. The small 
bubbles near + 30 m/s and - 30 m/s in range cells 65 
and above are indicative of returns at ranges where 
the return is weak with very little specfular structure. 
Investigation is continuing to use these methods 
coupled with an expert mode classification algorithm 
to detect hazardous windshear. 

Ref: M.W. Kunkel, "Spectrum Modal Analysis for the 
Detection of Low-Altitude Windshear with Airborne 
Doppler Radar," Radar Systems Lab TR-15, ECE Dept., 
Clemson University, Feb. 1992. 
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NOTES 


This is a ground track plot of one of several clutter 
flights over runway 26L at Denver Stapleton. 
Generally these flights were at altitudes of 1400 
to 1500 feet at groundspeeds near 200 knots. Various 
antenna elevations were used. Results here include 
data from horizontal elevation ( 0 degrees) and -1 and 
-3 degrees (below horizontal). The position plot is 
relative to the location of the TDWR (0 0) 




Ground Clutter Doppler Spectrum (den rwy 26 Flyover) 
1440 ft alt., AZ=0, EL=-1 
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NOTES 

A typical clutter FFT spectrum is shown which includes 
"discrete clutter" modes at about -18 m/s and +25 m/s 
in addition to the main lobe clutter near zero Doppler. 
These modes away from zero are due to returns from 
interstate highway 1-70 which passes under runway 35. 
Also shown is the impulse response of a second order 
adaptive clutter rejection filter with the I &Q sequence 
used for the illustrated spectrum also used as a training 
sequence for the adaptive filter. Notice that the adaptive 
filter places not ches at each of the strong clutter modes 
even with only a second order filter. The filter is an 
adaptive noise canceller using the LMS algorithm and 
is shown in block diagram form in the next figure. 
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Adaptive Noise Canceling with LMS Filter 


LMS Algorithm: *,(/» + 1) = w,(n) +^e(n)yj(n - /), 0 si 
where e(n) = x(n) + v,(n) - g(n ) = x(n) + _ w. r v 2jl . 


Primary 



Tapped-delay-lioe model 
Adaptive noisc-cancclcr configuration summary. 
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NOTES 


An adaptive noise canceller can be used to optimally 
reject clutter from a radar return if a reference clutter 
sequence is available and that sequence is highly 
correlated with the ground clutter po rtion of the 
primary return and uncorrelated with the weather 
portion of the primary return. As shown earlier a 
low order tapped delay type adaptive filter may be 
very effective and is being investigated as an 
alternative for very low signal to clutter ratio situations. 
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Analysis of Modelled Microbursts in Real Clutter 
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NOTES 

To investigate very low signal to clutter ratio situations 
the NASA simulation model is being used in conjunction 
with actual recorded ground clutter returns. The 
simulation model is set up with a microburst windfield 
from a previously observed microburst. It can be 
placed at any location for which clutter data have 
been recorded. Simulated I & Q data are then simply 
added to the archived clutter I & Q data after proper 
scaling and used for analysis of signal processing 
algorithms. The weather signal to clutter ratio can 
be controlled as desired by the user. 
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Adaptive Noise Canceling Clutter Rejection Results 

post clutter rejection filter pulse-pair mean estimate error 
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NOTES 

The LMS based adaptive noise canceller described 
earlier was used to investigate detection of the 
July 11 1988 Denver microburst in the presence 
of actual radar ground clutter data recorded in 1991. 
The weather return levels were scaled to maintain 
a constant average signal to clutter ratio of -20 dB 
in each range cell. The pulse pair estimate of Doppler 
spectrum mean was then computed after filtering 
with an adaptive noise canceller using the recorded 
and time delayed clutter data as a reference input. 
Results were compared to similar processing after 
filtering with a fixed 1.5 m/s notch Butterworth filter. 
The adaptive filter is much better where discrete 
clutter interference is present (near range cells 50 
and 64). This should improve hazard factor estimates. 




DEN rwy 26 Flyover, el=o, 3755 prf, 1-1 4km 
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NOTES 


This and the next 8 slides analyze the effect of out-of -range 
returns in the Denver data when the aircraft was flying 
over runway 26 headed west with the radar antenna 
oriented toward the mountains.Three views show the 
radar Doppler spectra for range cells 6-96 at the 3755 
prf. Data at 1877.5 prf were also recorded and analyzed 
to show what is in range cells 6-96 without second 
time around returns and what is in the extended range 
cells that aliases into the closer range cells when 
operating at the higher prf. Three different antenna 
elevation angles are shown. Out-of-range interference 
is significant at the horizontal elevation (EL=0) and is 
reduced to a negligible level at EL=-3 degrees. Some 
spectrum aliasing is also noted in the reduced prf 
plots since the Doppler range is halved to -15.+15 m/s. 
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DEN rwy 26 Flyover, el=-i, 3755 prf, i-i4km 
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DEN rwy 26 Flyover, el=-i, i 877.5 prf, 42-56km 
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DEN rwy 26 Flyover, el=-3, i 877.5 prf, i-i4km 

dn3clm3 Fromm 210 re 6-96 16 .ompl.. PRF/2 
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Mean Goundspeed Error from Radar Return - DEN 
dn2c6as1.m6, EL=0, all times with corrected AZ * 0, rc 6-96 
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NOTES 


One of the Denver runway 26 overrflights is analyzed 
to determine any difference from zero of the Doppler 
location of the ground clutter mainlobe. This difference 
is interpreted as a difference between the groundspeed 
used to determine Doppler zero in the radar 
demodulation and the groundspeed as measured by 
the radar. In each frame of 128 radar pulses the 
spectrum peaks in all range cells have been averaged 
to get a "frame" mean. These are then plotted for 
each frame when the corrected antenna azimuth 

was 0 degrees. A slight bias of aproximatedly 0.055 m/s 
is noted. 



St. Dev. Goundspeed Error from Radar Return - DEN 
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NOTES 


The situation described in the previous slide has 
been analyzed to estimate the standard deviation 
in the "groundspeed error" as determined from 
the radar measurement. 
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Groundspeed Correction Bias vs. Antenna Azimuth 
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NOTES 

During each flight experiment the airborne radar 

a il! en I!? was in a scanning mode, typically coverina 
-30, +30 degrees azimuth relative to the A/C longitudinal 
axis. A small bias error in the ground speed that 
varied linearly as a function of azimuth angle was 
noted in yirtually all data sets. This error remained 
after routine correction for the geometric variation 
0 ^°PP* er ss a function of azimuth angle. In 
addition, the slope of this bias varied depending on 
the antenna scan direction. This figure illustrates 
a counterclockwise scan from a DEN rwy 26 flyover 
FFT spectra of data records from 61 frames at 1 
degree increments across the scan were thresholded 
so that the central ground clutter peak is represented 
as a cluster of points showing how the Doppler spectrum 
peak vanes with azimuth angle. All data are from the 

range cell where the antenna boresight intersects 
the ground. 




Additional On-going Research Work 


1. Real time algorithm development 


2. Adaptive d utter rejection filtering 


3. Modelling of weather radar returns 

- autoregressive modelling 

- linear models based on fractals 

- non-linear models based on chaotic systems theory 


4. Optimal Bayesian methods for Hazard Factor estimation 

- use Doppler spectra as windspeed probability density 

- transform to F-factor probability density 

- estimate ‘most probable' F-factor map in protection volume 
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